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ABSTRACT 


r 

The  vertical  plane  motions  of  SWATH  ships  are  theoreti¬ 
cally  modeled.  Strip  theory  is  used  to  evaluate  hydrodynamic 
forces.  Contributions  due  to  body  lift,  cross  flow  drag,  and 
fin  lift  dominate  the  damping  coefficients.  Consequently, 
their  accurate  modeling  is  vital  to  the  accuracy  of  motion 
predictions.  Semiempir ical  methods  developed  for  evaluating 
these  components  are  described.  Data  for  oscillating  two- 
dimensional  cylinders,  flat  plates,  and  pairs  of  fins  as 
well  as  semiempirical  expressions  for  submarine  hydrodynamic 
coefficients  have  been  utilized  in  this  development.  Correla¬ 
tion  between  predicted  and  experimental  results  are  presented 
for  hydrodynamic  coefficients,  exciting  force  and  moment,  and 
responses  to  regular  waves.  The  expressions  developed  result 
in  correlation  which  is  good  and  which  is  notably  better  than 
results  reported  previously.. 

ADMINISTRATIVE  INFORMATION 

This  work  was  funded  under  the  Ships,  Subs  and  Boats  Program  Task  Area  SF 
4il-350-200,  N62345.  The  funding  was  administered  by  the  Exploratory  Development 
Programs  Office,  Code  1506,  Ship  Performance  Department,  David  Taylor  Naval  Ship 
Research  and  Development  Center  (DTNSRDC) . 


INTRODUCTION 

The /Small  Waterplane  Area  Twin  Hull  (SWATH)  ship  is  composed  of  two  hulls, 
each  of  which  has  one  or  two  surface  piercing  struts  connecting  them  to  the  above 
waterline  deck.  Typically  the  lower  hull  is  composed  of  circular  or  elliptical 
cr  s-sections.  Some  hulls  are  submarine-like  in  shape  and  others  are  composed 
of  a  series  of  cylinders  and  conic  frustums. 

The  motions  of  SWATH  ships  are  greatly  determined  by  the  ship's  unique  geome 
try.  Since  a  large  portion  of  the  ship's  buoyancy  is  located  in  the  lower  hulls, 
the  wave  exciting  forces  are  relatively  small.  The  waterplane  area  (Aw)  and 
longitudinal  metacentric  height  (GM^)  are  small  in  comparison  with  those  of  con¬ 
ventional  displacement  ships.  Since  the  heave  and  pitch  natural  periods  are 
inversely  proportional  to  the  square  root  of  Aw  and  GM^,  respectively,  relatively 
long  natural  periods  result.  Low  responses  for  operation  at  moderate  speeds  in 
seaways  occur  in  part  because  most  of  the  energy  of  a  seaway  typically  occurs  at 
short  wave  periods.  Thus,  the  small  waterplane  area  can  result  in  a  highly 
seaworthy  ship. 


Two -dimensionni  theory  has  been  used  to  predict  the  motions  in  a  seaway  of 
1*  2 

conventional  ships  and  SWATH  ships.  Two-dimensional  theory  assumes  that  there 
is  no  longitudinal  hydrodynamic  interaction  so  that  hydrodynamic  forces  can  be 
evaluated  by  integrating  the  hydrodynamic  contributions  of  two-dimensional  sections 
along  the  ship's  length. 

For  conventional  displacement  ships  and  for  catamarans,  as  well,  vertical 

plane  motions  can  be  predicted  accurately  using  potential  flow  theory.  However, 

2 

Lee  recognized  that  for  SWATH  ships  viscous  contributions  to  the  hydrodynamic 

damping  coefficients  are  important.  They  can  be  dominant,  making  their  accurate 

2 

modeling  important.  In  a  theoretical  development  Lee  introduced  contributions 

3 

due  to  lift  and  cross  flow  drag  of  the  body  and  stabilizing  fins.  Hong  introduced 
pitch  due  to  surge  and  demonstrated  its  importance  in  modeling  low  speed  motions. 
This  approach  was  generally  successful  in  predicting  the  vertical  plane  responses 
of  SWATH  ships.  However,  discrepancies  between  predicted  and  experimental  magni¬ 
tudes  and  uncertainty  ever  the  appropriate  values  of  cross  flow  drag  and  lift 
coefficients  motivated  the  present  study.  The  goal  of  this  investigation  is  to 
improve  the  quality  of  predictions  and  to  predict  vertical  plane  responses  of  a 
SWATH  ship  to  waves  given  only  the  ship  geometry,  the  location  of  the  center  of 
gravity,  and  the  longitudinal  radius  of  gyration  (gyradius). 


EQUATIONS  OF  MOTION 

For  the  purpose  of  this  derivation,  a  SWATH  ship  is  assumed  to  be  moving  at 
a  constant  speed  at  a  fixed  angle  relative  to  a  sinusoidal  wave  train  in  infinitely 
deep  water.  The  wave  amplitude  is  assumed  to  be  small  so  that  the  rigid  body 
motions  can  be  described  using  a  linear  model.  The  ship  is  defined  in  a  right- 
handed  coordinate  system  having  its  origin  at  the  mean  waterline  at  the  ship's 
longitudinal  center  of  gravity  and  centerline.  The  z-ordinate  is  positive  upward. 

The  equations  of  motion  of  the  ship  for  the  vertical  plane  are: 


MCi 
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where  M  is  the  mass  of  the  displaced  volume  and  is  the  pitch  mass  moment  of 
inertia.  A..,  B . . ,  and  C..  are  the  added  mass,  damping,  and  restoring  coefficients 

*.  L  J  J  J  A-  I 

in  the  i  mode  due  to  a  sinusoidal  motion  of  unit  amplitude  in  the  j  mode.  The 
subscript  1  denotes  surge,  3  heave,  and  5  pitch.  is  the  complex  exciting 

force  or  moment  and  gj  is  the  wave  frequency  of  encounter. 

In  evaluating  the  coefficients,  forces,  and  moments,  strip  theory  is  employed. 
This  is  a  reasonable  assumption  for  the  SWATH  ship  with  its  slender  and  gradually 
changing  geometry  over  its  length. 


HYDRODYNAMIC  COEFFICIENTS  AND  EXCITING  FORCES  AND  MOMENT 

The  hydrodynamic  coefficients  and  exciting  forces  and  moment  are  composed  of 

2 

potential  flow,  cross  flow  drag,  and  lift  terms.  Lee's  development  included 
cross  flow  drag  and  lift  terms  for  the  body  and  the  stabilizing  fins.  These 
components  affect  the  damping  and  restoring  coefficients  and  the  exciting  forces 
and  consequently  the  motions  at  all  speeds  since  the  cross  flow  drag  terms  are 
independent  of  speed  and  the  lift  terms  are  proportional  to  some  power  of  speed. 

In  this  investigation,  new  expressions  evolved  for  the  cross  flow  drag  coeffi¬ 
cients,  the  body  lift  terms,  and  the  fin  cross  flow  drag  and  lift  curve  slope 
coefficients.  Final  expressions  for  the  hydrodynamic  coefficients  and  the  exciting 
forces  and  moments  are  given  in  Appendix  A.  Details  of  the  development  of  these 
expressions  are  given  in  Appendices  B  and  C.  It  is  useful  to  briefly  summarize 
the  results. 

Currently,  the  potential  flow  components  of  the  added  mass  and  damping  can  be 

4 

evaluated  utilizing  either  the  Frank  Close  Fit  Technique  or  the  Dalzell  Approxima¬ 
tion  Technique.*  When  the  Dalzell  Approximation  Technique  is  utilized,  as  it  is 
in  the  results  in  this  report,  the  distribution  of  the  potential  on  the  two- 
dimensional  section  is  unavailable  and  evaluation  of  the  exciting  force  and 


*As  described  by  Dalzell  in  Stevens  Institute  of  Technology  reports  with 
limited  distribution. 
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moment  requires  further  approximation.  In  Appendix  B  expressions  for  the  surge 
and  heave  exciting  forces  and  the  pitch  exciting  moment  are  developed.  This 
development  includes  an  expression  for  an  approximate  depth  used  in  the  heave 
exciting  force  and  pitch  exciting  moment  and  an  expression  for  the  center  of  surge 
which  facilitates  inclusion  of  surge  in  the  exciting  moment. 

In  Appendix  C  a  development  of  the  cross  flow  drag  and  lift  terms  and  the 

corresponding  coefficients  is  given.  The  cross  flow  drag  coefficients  for  the 

body  are  evaluated  using  experimental  values  for  oscillating  circular  cylinders.^ 

A  factor  to  reflect  the  effect  of  a  strut  on  the  cross  flow  drag  coefficient  is 

developed  from  oscillation  data  for  two-dimensional  SWATH  sections.*  The  cross 

flow  drag  coefficient  for  the  stabilizing  fins  is  evaluated  using  experimental 

data  for  oscillating  plates.^  Semiempirical  expressions  for  the  vertical  plane 

hydrodynamic  coefficients  of  submarines**  serve  as  a  basis  for  the  development  of 

the  SWATH  body  lift  components.  The  lift  curve  slope  of  the  fins^  is  given  and  a 

correction  is  made  for  the  frequency  dependent  interference  effect  of  the  forward 

8  9 

fin  on  the  lift  of  the  aft  fin.  ’  For  appropriate  configurations,  an  additional 
correction  is  made  for  the  effect  of  the  hull  wake  on  the  lift  of  the  aft  fin.10 

COMPARISON  BETWEEN  RESULTS  FROM  EXPERIMENT  AND  PREDICTION 

Oscillation,  excitation,  and  regular  wave  data  are  available  for  SWATH  con¬ 
figurations  denoted  6A,  6B,  6C,  6D,  and  SSP  KAIMALINO.  These  results  were  used 
to  guide  the  development  of  the  expressions  in  this  report.  The  four  configura¬ 
tions  in  the  SWATH  6  series  employ  the  same  lower  hull  which  is  a  body  of  revolu¬ 
tion.  The  struts  are  designed  so  that  the  GM^  differs  for  each  strut  design. 
Configurations  6A  and  6B  are  single  strut  designs,  whereas  6C  and  6D  are  twin 
strut  designs.  Particulars  of  the  configurations  and  of  the  lifting  surfaces  are 
given  in  Tables  1  and  2. 

Correlation  between  experiment  and  prediction  are  given  in  Figures  1  through 
8.  Nondimensionalization  factors  for  added  mass,  damping,  exciting  force,  and 

* 

As  described  by  Stahl  in  a  DTNSRDC  report  with  limited  distribution. 

** 

As  described  by  Dempsey  in  a  DTNSRDC  report  of  higher  classification. 


4 


exciting  moment  are  given  in  Table  3.  Open  symbols  are  used  for  experimental  re¬ 
sults  and  solid  symbols  are  used  for  predicted  results  using  the  expressions  given 

2 

in  Appendix  A.  Predicted  results  using  Lee's  expressions  are  given  in  solid  lines 

3 

on  some  of  the  added  mass  and  damping  and  all  of  the  regular  wave  figures.  Hong's 
modifications  are  included  in  the  regular  wave  results. 

In  solving  the  cross  flow  drag  contributions,  it  is  necessary  to  know  the 
motion  of  the  model  relative  to  the  incoming  wave.  For  regular  wave  motion  pre¬ 
dictions,  this  component  must  be  solved  iteratively,  until  the  responses  of  the 
craft  converge;  that  is,  the  difference  between  the  estimated  and  computed  re¬ 
sponses  diminish  to  acceptable  values.  However,  in  the  case  of  forced  oscillation 
experiments,  the  motion  of  the  body  is  known  and  there  is  no  incident  wave.  Con¬ 
versely,  for  wave  exciting  experiments,  the  body  is  held  rigid  and  the  motion  of 
the  wave  is  known.  Therefore,  calculation  of  these  components  is  straightforward. 

ADDED  MASS  AND  DAMPING 

Two  sets  of  data  for  heave  and  pitch  forced  oscillation  tests  are  available 
for  the  SWATH  6A.  Results  from  a  1:31.2  scale  bare  hull  model  for  speeds  corre¬ 
sponding  to  full-scale  speeds  of  10,  20,  and  35  knots11  are  presented  in  Figure  1 
along  with  predicted  results.  Results  from  a  1:22.5  scale  model  with  and  without 

stabilizing  fins  for  speeds  corresponding  to  full-scale  speeds  of  0,  20,  and  28 
12 

knots  are  presented  in  Figure  2.  Included  are  predicted  results  based  on  the 
expressions  in  Appendix  A  and  Lee's  predicted  results  which  were  presented  in 
Reference  12.  These  latter  predicted  results  do  not  include  the  cross  flow  drag 
contributions.  Inclusion  of  these  terms  would  increase  the  magnitudes  of  the 
damping  terms,  most  significantly  at  zero  speed. 

It  is  useful  to  compare  experimental  results  given  in  Figure  la  (Fn  =  0.384) 
with  those  in  Figure  2e  (bare  hull).  It  is  expected  that  these  results  should  be 
close  in  value  since  they  are  for  identical  conditions.  Only  the  model  scales 
differ.  However,  the  results  for  A^j  given  in  Figure  la  are  smaller  than  those 
given  in  Figure  2e.  The  difference  in  results  for  A^  is  particularly  important 
at  higher  speeds  as  can  be  seen  in  Figures  la  and  2i.  The  predicted  results  agree 
well  with  the  measured  results  in  Figure  la;  however,  no  explanation  of  the 
difference  in  experimental  results  is  proposed  in  this  presentation. 


WAVE  EXCITING  HEAVE  FORCE  AND  PITCH  MOMENT 
13 

Fein  and  Stahl  carried  out  experiments  to  measure  the  surge  and  heave 
wave  exciting  forces  and  the  pitch  exciting  moment.  They  investigated  five 
speeds  in  head  and  following  seas  for  a  1:22.5  scale  model  of  the  SWATH  6D  and 
five  speeds  in  head  waves  for  a  1:7.8  scale  model  of  the  SSP  KAIMALINO.  The 
data  given  in  this  report  are  presented  in  a  different  format  from  that  of 
Reference  13.  To  elucidate  the  data  in  following  seas,  all  data  have  been 
presented  as  a  function  of  wavelength  to  ship  length,  rather  than  encounter 
frequency.  Since  the  theory  is  developed  with  the  pitch  moment  about  the  LCG 
at  the  mean  waterline,  the  measured  pitch  exciting  moment  and  surge  exciting 
force  were  used  to  transform  the  moment  to  be  about  the  LCG  at  the  mean 
waterline,  so  that  the  predicted  and  experimental  results  were  comparable. 

These  results  are  given  in  Figures  3  and  4. 

RESPONSES  TO  REGULAR  WAVES 
14  15 

Kallio  *  carried  out  regular  wave  experiments  for  the  SWATH  6  series. 

Heave,  pitch,  and  relative  bow  motion  responses  as  a  function  of  wavelength  to 

ship  length  are  given  for  five  relative  wave  headings  for  the  6A,  6B,  and  6C  and 

for  head  and  following  waves  for  the  6D.  Results  are  given  in  Figures  5  through 

8.  Note  that  two  sets  of  predicted  results  are  given  for  all  conditions.  One 

2 

set  results  from  the  development  in  this  report  and  one  results  from  Lee's  work 
3 

with  Hong's  modifications  included. 

DISCUSSION 

Initial  work  by  Lee  indicated  that  correlation  between  experimental  and  pre- 

3 

dieted  results  at  zero  speed  for  long  wavelengths  was  not  satisfactory.  Hong's 
results  demonstrated  the  importance  of  introducing  the  effects  of  surge  and  of 
using  the  proper  wave  amplitude  in  evaluating  the  nonlinear  terms.  The  improved 
correlation  that  results  from  the  incorporation  of  the  expressions  developed  in 
this  report  is  evidently  due  to  the  method  of  evaluating  Cp.  Previously,  it  had 
been  assumed  to  be  constant,  whereas,  for  two-dimensional  sections  (or  stabilizing 
fins)  it  is  here  considered  to  be  a  function  of  the  major  and  minor  axis  of  the 
lower  hull  and  the  strut  thickness  (or  fin  span),  the  wave  frequency  and  the  ampli¬ 
tude  of  the  relative  vertical  velocity  at  each  two-dimensional  section  (or  fin). 

High  speed  responses  for  all  configurations  utilizing  the  expressions  derived 
in  this  report  result  in  excellent  correlation  with  experimental  data.  The 


responses  predicted  through  the  present  methodology  are  notably  closer  to  the 
experimentally  determined  responses  than  results  found  with  any  of  the  previous 
methods. 

Responses  of  SWATH  configurations  traveling  at  high  speed  in  following  waves 

has  been  a  topic  of  interest.  One  problem  in  correlation  for  this  condition  is 

14 

that  it  is  difficult  experimentally.  As  noted  by  Kallio,  in  quartering  and 

following  waves  there  was  considerable  surge  and  consequently  the  model  safety 

2 

restraint  lines  became  taut.  Lee  reported  extremely  large  predicted  heave 

responses  for  the  6A  at  the  wavelength  corresponding  to  zero  encounter  frequency.* 

Since  the  potential  flow  two-dimensional  approach  is  not  valid  at  small  encounter 

frequencies,  theoretical  work  was  undertaken  to  overcome  this  limitation. 

Hong applied  to  SWATH  ships  unified  slender  body  theory  developed  by  Newman 
18 

and  Sclavounos.  Results  from  Hong's  implementation  did  not  improve  correlation 

and  did  not  remove  the  spike  in  the  6A  heave  predictions.  However,  the  results 

developed  here  which  focused  on  the  viscous  components  but  retained  the  two- 

2 

dimensional  potential  flow  approach  utilized  by  Lee,  do  not  include  the  spiked 
response.  Although  pitch  is  overpredicted  for  the  6A  and  6B,  these  results  show 
generally  good  correlation  and  support  the  hypothesis  that  the  aberrant  predictions 
for  the  6A  are  not  related  to  the  two-dimensional  potential  flow  theory. 

Since  the  predicted  heave  response  spike  occurs  near  zero  encounter  frequency, 
it  has  been  assumed  that  the  problems  with  the  predictions  were  due  to  two- 
dimensional  theory.  However,  Figure  5e  suggests  an  alternative  explanation.  The 
large  heave  response  reported  by  Lee  occurs  near  the  wavelength  which  corresponds 
to  zero  encounter  frequency;  however,  this  also  occurs  in  the  region  where  the 
pitch  response  peaks.  Since  heave  and  pitch  are  coupled,  errors  in  modeling  one 


•ft 

The  regular  wave  predictions  which  are  attributed  to  Lee  include  a  modifica¬ 
tion  which  was  an  attempt  to  remove  the  spiked  behavior.  In  the  modification  when 
w  is  less  than  0.07  the  potential  flow  added  mass  and  damping  coefficients  for  each 
section  have  been  assumed  to  be  equal  to  those  for  w  *  0.07.  Generally  this  will 
alter  predictions  only  at  high  speeds  in  following  or  stem  seas.  This  approxima¬ 
tion  merely  suppressed  the  response.  (See  Figure  5e.) 
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r 


response  will  affect  the  other.  That  is,  the  spike  in  the  predicted  heave  for  the 
6A  may  be  due  to  Inadequacies  in  the  pitch  predictions. 

Analysis  of  the  relative  Importance  of  various  terms  in  the  pitch  equation  of 
motion  indicates  that  the  term  may  be  dominant.  For  simplicity,  consider  the 
uncoupled  pitch  equation  of  motion: 

d5  +  A55)£5  +  B55£5  +  C55?5  “  F5e 

Utilizing  the  relationship  £,.  *  (C^R  +  i£^j)e”:I'a)t,  this  becomes 

[-(o2  (I5  +  A55)  +  C55]£5  +  B55£5  -  F5e_ia)t 

Two-dimensional  theory  is  not  valid  at  small  encounter  frequencies  and  Arc  becomes 

2  33 

very  large  in  this  region.  However,  the  presence  of  o>  results  in  a  very  fortunate 

2 

situation  for  small  w.  That  is,  when  03  approaches  zero,  (i)  A^^  will  be  small. 

The  important  term  is  C^.  Whereas  A^,.  and  are  calculated  using  two-dimen¬ 
sional  theory,  C,^  is  not.  is  composed  of  a  fin  lift  term,  a  body  lift  term, 

and  a  term  which  is  essentially  GM^  Since  the  fin  lift  term  will  be  approximately 
equal  for  all  configurations  in  the  SWATH  6  series,  it  can  be  neglected  in  this 
discussion.  In  addition,  as  configured  in  this  report,  the  body  lift  term  is 
dependent  on  mass  and  particulars  of  the  lower  hull  and  will  be  equal  for  the  6A, 
6B,  and  6C  configurations.  However,  GM^  increases  significantly  from  the  first 
model  in  the  SWATH  6  series  to  the  last  with  the  6A  having  the  smallest  and  the 
6D  the  largest  GM^.  For  small  GM^  the  lift  terms,  and  the  body  lift  term  in 
particular,  will  be  relatively  more  important  than  for  a  large  GM^  configuration. 
This  argument  is  consistent  with  the  correlation  which  indicates  that  the  pre¬ 
dicted  heave  spike  in  earlier  work  occurred  for  the  6A  only.  Consequently,  good 
correlation  between  experimental  and  predicted  results  for  configurations  with 
small  GM^'s  will  be  strongly  dependent  on  accurate  modeling  of  C^. 

The  body  lift  component  used  for  the  development  is  based  on  experimental 
work  on  the  6A.  This  is  an  unfortunately  slim  data  base.  It  is  expected  that 
this  term  is  important  for  other  low  GM^  configurations.  Better  correlation  for 
pitch  for  the  6A  and  6B  probably  would  result  from  better  modeling  of 

! 
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CONCLUSIONS  AND  RECOMMENDATIONS 


2 

1.  Work  by  Lee  has  been  used  as  the  basis  for  modeling  the  vertical  plane 
motion  of  SWATH  ships.  The  general  form  of  his  work  has  been  retained  in  the 
present  modeling;  however,  the  effort  reported  here  has  focused  on  the  viscous 
terms.  Alternate  semiempirical  expressions  for  body  lift  terms,  cross  flow  drag 
coefficients  for  the  body  and  the  fins,  and  lift  curve  slope  coefficients  for  the 
fins  have  been  developed,  and  improved  correlation  with  experiment  is  the  result. 

2.  Correlation  of  added  mass  and  damping  coefficients  and  the  exciting  force 
and  moment  are  generally  good.  Zero  speed  response  correlation  is  comparable  to 
or  better  than  results  based  on  Hong's  modifications  to  Lee's  work.  High  speed 
correlation  is  good  and  is  notably  better  than  previous  results  shown  by  Hong.^ 
Following  sea  results  no  longer  display  the  aberrant  behavior  which  occurred  in 
Lee's^  and  Hong's^  SWATH  6A  results. 

3.  The  two-dimensional  potential  flow  theory  is  certainly  adequate  for  the 
prediction  of  ship  motions  of  SWATH  ships  similar  to  the  6  series. 

4.  Further  experimental  work  would  facilitate  refinement  of  the  viscous 
expressions  developed  here  and  would  expand  the  regions  of  confidence.  The  follow¬ 
ing  investigations  are  recommended: 

a.  Experimental  investigation  of  C^  for  SWATH  sections  and  for 
circular  cylinders  at  very  low  Kc  would  be  useful. 

b.  Experimental  investigation  of  the  effect  of  fin-fin  interference  for 
additional  configurations,  including  ones  where  the  aft  fin  is  larger 
than  the  forward  fin  would  be  useful. 

c.  Experiments  of  true  f ree-to-surge  conditions  using  radio-controlled 
models  in  following  waves  would  aid  in  the  assessment  of  prediction 
techniques. 

d.  Experimental  investigation  of  C55  for  existing  models,  including  the 
6B,  6C,  and  6D,  would  make  it  possible  to  define  the  body  lift  component 
of  more  precisely  and  to  improve  the  reliability  of  prediction  of 
responses  in  following  seas  at  high  speeds. 
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TABLE  1  -  FULL-SCALE  PARTICULARS  OF  SWATH  CONFIGURATIONS 


Particular 

6A1 

6B1 

6C1 

6D2 

3 

SSP 

Length  overall,  m 

73.15 

73.15 

73.15 

73.15 

26.40 

Distance  between  centerlines,  m 

22.86 

22.86 

22.86 

26.80 

12.19 

Draft,  m 

8.13 

8.13 

8.13 

8.13 

4.66 

Displacement,  metric  ton 

2,946 

2,946 

2,946 

2,946 

2,946 

Longitudinal  CG,  aft  of  nose,  m 

35.45 

35.14 

34.72 

36.10 

13.46 

Vertical  center  of  gravity  (KG),  m 

10.36 

10.36 

10.36 

9.00 

4.28 

Longitudinal  metacentric  height,  m 

6.10 

11.60 

13.70 

26.40 

5.03 

lAs  given  in  Reference  14. 
*As  given  in  Reference  15. 


As  given  in  Reference  13. 
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TABLE  2  -  PARTICULARS  OF  STABILIZING  FINS 


Configuration 

6A1 

6C1 

6D2 

3  4 

SSP  ’ 

Forward  Fin,  Each 

Chord,  m 

2.59 

2.16 

2.16 

2.59 

1.95 

Span,  m 

3.11 

2.59 

2.59 

3.11 

1.83 

Location,"*  m 

17.15 

17.15 

17.15 

17.15 

2.82 

Aft  Fin 

Chord,  m 

4.48 

3.73 

3.73 

4.48 

2.38 

Span,  m 

5.36 

4.48 

4.48 

5.36 

10.55 

Location,"*  m 

62.24 

62.24 

62.24 

62.24 

20.59 

^As  given  in  Reference  14. 

2 

As  given  in  Reference  15. 

3 

As  given  in  Reference  13. 

4 

Aft  fin  spans  between  the  hulls. 

^Longitudinal  distance  between  lower  hull  nose  and  fin  quarter  chord 
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TABLE  3  -  NONDIMENSIONALIZATION  FACTORS  FOR  ADDED  MASS, 
DAMPING,  WAVE  EXCITING  FORCE,  AND  WAVE  EXCITING  MOMENT 


Variable 

Nond imen  s iona 1 iza  t ion 

Factor 

A33 

M 

a35,  a53 

ML 

A55 

ML2 

B33 

M(g/L)1/2 

b35,  b53 

M(gL) 1/2 

B55 

ML(gL) 1/2 

F<«> 

MgA/L 

F<*> 

MgA 

O) 

(g/L)1/2 

where  A 

=  Wave  mplitude 

g 

=  Acceleration  due  to  gravity 

L 

=  Overall  ship  length 

M 

=  Mass  of  Displaced  volume 
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Figure  1  -  Comparison  between  Experiment  and  Prediction  for  Added 
Mass  and  Damping  Coefficients  of  the  SWATH  6A  (Bare  Hull) 

for  Various  Speeds 


Figure  1  (Continued) 
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Figure  1  (Continued) 
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Figure  lc  -  Damping  Coefficients  and 
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Figure  2  (Continued) 
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Figure  2  (Continued) 


Figure  2e  -  Added  Mass  Coefficients  A, 
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Figure  2  (Continued) 


Figure  2g  -  Damping  Coefficients  B 


for  Fn 


Figure  2  (Continued) 
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Figure  2  (Continued) 
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Figure  2  (Continued) 


Figure  2k  -  Damping  Coefficients  and  for  Fn  = 


0.52 


Figure  2  (Continued) 


Figure  21  -  Damping  Coefficients  Bi-  and  B'  for  Fn  =  0.53 
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Figure  3  -  Comparison  between  Experiment  and  Prediction  of  Heave 
Exciting  Force  and  Pitch  Exciting  Moment 
for  the  SWATH  6D 
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Figure  3a  -  Heave  Force  for  Fn  =  0.0  in  Head  Waves 
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Figure  3  (Continued) 
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Figure  3b  -  Pitch  Moment  for  Fn  *  0,0  in  Head  Waves 


Figure  3  (Continued) 
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Figure  3c  -  Heave  Force  for  Fn  =  0.077  and  0.192  in  Head  Waves 
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Figure  3  (Continued) 
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Figure  3d  -  Pitch  Moment  at  Fn  =  0.077  and  0.192  in  Head  Waves 


Figure  3  (Continued) 
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Figure  3e  -  Heave  Force  at  Fn  *  0.384  and  0.538  in  Head  Waves 


37 


►  ■  ♦  •  • 
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Figure  3  (Continued) 


Figure  3g  -  Heave  Force  for  Fn  =  0.0  in  Following  Waves 
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Figure  3  (Continued) 
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Figure  3  (Continued) 


Figure  3k  -  Heave  Force  at  Fn  =  0.384  and  0.538  in  Following  Waves 
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Figure  3  (Continued) 


Figure  3d  -  Pitch  Moment  at  Fn  =  0.384  and  0.538  in  Following  Waves 
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Figure  7  -  Comparison  between  Experiment  and  Prediction  of 
Regular  Wave  Transfer  Functions  for  the  SWATH  6C 
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APPENDIX  A 


FREQUENCY  DEPENDENT  HYDRODYNAMIC  COEFFICIENTS,  FORCES  AND  MOMENTS 
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where  integrations  are  along  ship  length  with  the  origin  at  the  LCG.  Contributions 
from  both  hulls  are  included  in  the  integrations. 


AR  = 


a  +  2s 


A  =  Projected  area  of  n^  stabilizing  fin 


Aw  =  Waterplane  area 


a  = 
o 


-L2(2.439((k2  -  kj)m')4/3  -  m') 
f  a  dx 

Ja 


a^^  =  Heave  added  mass  coefficient  of  two-dimensional  section 

it  2 

a33  =  Added  mass  of  stabilizing  fins  P sc 

n 


b^^  =  Heave  damping  coefficient  of  two-dimensional  section 
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(see  Figure  9  for  C 


S ARP KAY A 


) 


SARPKAYA 


th 


Cp  for  n  stabilizing  fin  modeled  by  flat  plate  with  d  =  2s 
(see  Figure  10) 

1  8  7T  AR 

- 1— r - T-p.  C  (with  C  for  aft  fins  is  corrected 

1.8  +  (AR2  +  4 . 0) 1/2  n  .¥n  _  .  .  _  ..  . . 

according  to  data  m  Table  4) 


‘■°  -  -aVH'/r  -  °-1612] 


1/2 


0.6366  sin 


-1  0.4015 


(a  +  2s)/a 


if  for  aft  fin  >  2 


.1.0  otherwise 

/  Bt  V 

(l  ~  S  +  Bt  )  di 


a 

o 


2  -  a 

o 


8 

o 


2  -  6o 


=  Mass  of  displaced  volume  of  SWATH  ship 


M 


7T 

=  Mass  of  stabilizing  f in  =  p  (set) 

=  Total  number  of  fins 

■  T  ab 

=  Transverse  distance  between  fin's  centroid  and  the  ship's 
centerline 


t  =  Draft  of  strut 
s 
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Figure  9  -  Drag  and  Inertia  Coefficients  versus  Keulegan-Carpenter 
Number  for  Constant  Values  of  the  Frequency  Parameter 
(From  Reference  5) 
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Figure  10  -  The  Drag  Coefficient  of  Flat  Plate  (+) , 
Diamond  (0^  and  Circular  (Q>  Cylinders 
at  Low  KC  (From  Reference  6) 
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r  T 


TABLE  4  -  RATIO  OF  LIFT  ON  AFT  FIN  TO  LIFT  ON  FORWARD  FIN  FOR 
VARIOUS  FIN  SEPARATIONS  AND  OSCILLATION  FREQUENCY 
TO  SPEED  RATIOS  (FROM  REFERENCES  3  AND  9) 


K  \  GJS 

r\u 

0.00 

0.04 

0.08 

0.12 

0.16 

0.20 

10 

0.412 

0.544 

0.643 

0.824 

1.076 

1.221 

15 

0.462 

0.638 

0.846 

1.046 

1.180 

1.109 

20 

0.529 

0.732 

1.000 

1.151 

1.110 

0.971 

25 

0.614 

0.816 

1.099 

1.132 

1.011 

0.897 

20 

0.706 

0.853 

1.118 

1.006 

0.912 

0.853 

x  =  distance  between  leading  edges  of  fins 
s  =  span 

U  =  forward  speed 

cj  =  oscillation  frequency 


where 


APPENDIX  B 

POTENTIAL  FLOW  COMPONENTS 


In  Reference  2,  Lee  discusses  the  potential  flow  added  mass  and  damping 
coefficients  for  oscillating  twin  cylinders.  The  source  distribution  method  is 
applied.  This  technique  requires  a  substantial  amount  of  computer  time  and  con¬ 
sequently  limits  the  number  of  configurations  which  can  be  evaluated.  John 
* 

Dalzell  initiated  an  effort  to  approximate  the  potential  flow  added  mass  and 
damping  coefficients.  The  goal  of  his  work  was  to  develop  expressions  which 

would  approximate  the  potential  flow  coefficients  resulting  from  the  Frank  Close 

4  19 

Fit  Technique  which  is  utilized  in  References  1  and  2.  Lee  also  developed 

expressions  for  heave  added  mass  and  damping  and  extended  the  approximate  approach 

* 

to  the  transverse  plane.  Dalzell  then  modified  his  work  further  and  extended  it 
to  include  the  transverse  plane  and  elliptical  as  well  as  circular  lower  hull 
sections.  Dalzell' s  expressions  for  approximating  the  heave  added  mass  and 
damping  coefficients  for  two-dimensional  sections  have  been  used  in  results  in 
this  report.  His  work  correlates  well  with  results  from  the  Frank  Close  Fit  for 
a  variety  of  strut-lower  hull  configurations.  Consequently,  the  expressions 
developed  in  this  report  for  the  viscous  components  can  be  expected  to  result  in 
similar  values  when  used  in  conjunction  with  potential  flow  results  obtained  using 
Frank  Close  Fit  Technique  or  the  Dalzell  Approximation  Technique.  One  consequence 
of  using  the  Dalzell  Approximation  Technique  is  that  the  distribution  of  the 
potential  on  the  two-dimensional  section  is  unavailable,  necessitating  approxima¬ 
tions  in  the  exciting  forces  and  moment. 

The  potential  flow  exciting  force  and  moment  include  the  incident  wave 

potential  and  a  diffraction  potential.  Salvesen,  Tuck  and  Faltinsen^  show  that 

(e)  (e) 

the  exciting  force  and  moment  F^  and  F^  are 
F^  =  PA  J"  (f3(x)  +  hj (x) )dx 

F^6)  =  PA  J'  T  x (f  3 (x)  +  h3(x))  +  h3(x)l  dx 


As  described  by  Dalzell  in  Stevens  Institute  of  Technology  reports  with 
limited  distribution. 
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where  the  sectional  incident  wave  contribution,  which  is  often  referred  to  as  the 
Froude-Krylov  force,  is  defined  by 


f3(x) 


ge 


ikx  cosB 


/ 


n^e 


-iky 


sin  8  kz 
e 


dl 


and  the  sectional  diffraction  contribution  is  defined  bv 


h„  (x)  =  -0)  e 
J  o 


ikx 


cosB  J" 


(in^  +  n9sinB)e  sin^ekz,|^  dl 


where  p  is  the  mass  density  of  water,  A  is  the  wave  amplitude,  w  is  the  wave 
frequency,  (a)  is  the  wave  encounter  frequency,  k  is  the  wave  number,  8  is  the 
heading  of  the  ship  relative  to  the  incident  wave,  U  is  the  forward  speed  of  the 
ship,  n 2  and  n^  are  unit  normals,  is  the  velocity  potential,  and  dl  is  an 
element  of  an  arc  along  the  cross-section  C^.  (End-effect  terms  have  been  ex¬ 
cluded  here  since  they  are  unnecessary  for  SWATH  configurations.) 

20 

Korvin-Kroukovsky  made  the  assumption  that  an  average  depth  can  be  utilized 
lcz 

so  that  the  term  e  can  be  moved  outside  the  integral.  Salvesen,  Tuck  and 
Faltinsen  investigated  this  approach.  The  complex  potential  can  be  expanded  in 
a  Taylor  series.  If  exp(-iky  sinB)  is  expanded  through  order  k  and  symmetry  about 
the  y=0  axis  is  assumed  then 


ikx 

f3(x)  =  ge 


cos8  -kd^ 
e 


n^  dl 


ikx  cos 8  -kd^ 

=  ge  e  B 


where  B  is  the  local  waterline  beam  and  d3  is  an  appropriate  average  depth. 
Similarly,  h^(x)  becomes 


h3(x) 


P 


ikx  cosB  -kd 

woo  e  e  S(a_0  (u>) 

o  i  i 


b3 3  (ui)  ) 


(B.l) 


where  dg  is  an  average  depth  and  a^3  and  b33  are  the  sectional  added  mass  and 

damping  coefficients. 
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Newman  follows  a  different  approach.  In  his  evaluation  of  the  Froude- 
Krylov  force,  Newman  expands  the  exponential  under  the  integrand  giving 


f^(x)  =  cos&  j  n^(l  -  iky  sing  +  kz)dl 

^x 


ita  cosB[b(i()  .  ks(x)-| 


(B.  2 ) 


S(x)  is  the  local  sectional  area.  For  SWATH  sections  with  wall-sided  struts,  the 
argument  can  be  made  that  since  n^  is  zero  along  the  strut,  the  closed  curve  can 
be  approximated  by  the  ellipse.  This  approximation  which  provides  good  correla¬ 
tion  with  Frank  Close  Fit  results  has  been  utilized  here.  It  leads  to  S(x)  in¬ 
cluding  only  the  area  of  the  ellipse  for  SWATH  sections.  This  form  of  f^(x)  has 
an  advantage  over  the  previous  form  since  expansion  of  the  exponential  eliminates 
the  need  for  an  equivalent  depth. 

Application  of  Equations  (B.l)  and  (B.2)  for  h^(x)  and  f^(x)  makes  it 
possible  to  evaluate  the  exciting  heave  force  and  pitch  moment  if  the  sectional 
added  mass  and  damping  are  known  and  an  appropriate  average  depth  is  used. 
Solution  for  dg  using  numerical  results  from  the  Frank  Close  Fit  Technique  and 
from  Equation  (B.2)  for  SWATH  sections  indicates  that  dg  can  be  approximated  by 


Bt 


dg(x) 


1  + 


-jj-  ab  +  Btg 


where  B  is  the  beam,  t  is  the  draft  of  the  strut,  and  d.  is  the  distance  between 

s  1 

the  mean  waterline  and  the  center  of  the  ellipse. 

When  the  potentials  on  the  surface  of  the  body  are  unknown,  approximations 

to  the  surge  exciting  force  and  the  center  of  surge  force  are  necessary.  The 

(e) 

surge  exciting  force,  F^  ,  is  given  by 

_(e)  „  .  /*  ikx  cosS  -  iky  sin$  f  -kz 

F.  =  PgA  /  e  ^  j  e  n.  dx 

J0  •'c(x) 
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and  can  be  approximated  by 


ikx  co 
e 


dB  \ 
dx  ) 


dx 


Evaluation  of  the  expressions  in  Reference  3  for  various  SWATH  sections  using  the 
Frank  Close  Fit  Technique  led  to  an  approximation  for  d^,  the  center  of  surge 
force 


d  = 
e 


d 


1 


That  is,  the  exciting  moment  due  to  surge 
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APPENDIX  C 
VISCOUS  COMPONENTS 


2 

Lee  utilized  an  equation  described  in  Reference  22  to  model  the  viscous 
damping  of  SWATH  ships.  The  vertical  force  Fy  on  a  slender  moderately  inclined 
body  is  expressed  as 


Fv  “  T  PU\  sino|sina| (aQ|cota|  +  Cp) 

where  p  is  the  density  of  water,  U  is  the  ship  speed,  A^  is  the  projected  area  of 

the  body,  a  is  the  trim  angle,  aQ  is  the  viscous  lift  coefficient  and  Cp  is  the 

cross  flow  drag  coefficient.  For  air  ships  with  circular  cross  sections,  experi- 

22 

mental  work  indicates  that  a  is  about  0.07  and  Cn  is  between  0.4  and  0.7. 

2  o  D 

(Lee  suggests  using  0.07  and  0.5,  respectively.) 

Substituting  w/U  ■  sin a  and  assuming  small  a  this  becomes 

Fv  "  T  PAp(aou2°t  +  CDWIWI)  (C-1} 

Following  Lee,  w  is  the  relative  fluid  velocity  and  can  be  defined  as 

w  *  £3  -  xC5  +  y -Cv(x,±Sp(x),-d1(x),t) 
where  d^(x)  is  the  depth  of  submergence  of  the  maximum  breadth  of  a  section.  The 
relative  velocity  of  the  fluid  induced  by  the  incoming  wave,  Cv*  Is  given  by 

'  ,  .v  .  kz  +  ikx  cosB  -  iky  sinfi  -iwt 

Cv(x,y.z,t)  *  -iwoAe  7  e 

The  angle  of  incidence  of  flow  is  given  by 

O  -  S5  +  (C3  -  x£5  +  yC4  -  Cv(x,±Sp(x),d1(x)))/U 

The  equation  for  the  vertical  force  F  and  the  pitch  moment  M  are  assumed  to  be 


mmammmrn 


M  -  -J  xFy(x)dx 

where  Fv(x)  Is  given  by  Equation  (C.l). 

When  the  necessary  algebra  and  Integrations  along  the  body  have  been  per¬ 
formed,  and  the  equ 11 inear izat ion  approximation  has  been  applied  the  viscous 


contributions  become  : 

B33  "  P  ~W  CD  f  z10a(x)dx  +  T  aoU  f  a^x^dx  (C.2) 

B35  “  B53  "  “P  "ST  CD  /  xz^a(x)dx  -  -j-  aQU  xa(x)dx  (C.3) 

B55  "  P*5F  S)  J  x2z1Qa(x)dx  +  j-  a0u  J"  x2a(x)dx  (C.4) 

^  =  -j-  f  a(x)dx  (C.5) 


C5S  •  -  T  V2/ 


xa(x)dx 


(C.  6) 


*2  I  +  be  cosB)  ikSp  sin3( ,  -ikS^  sinB 


F;  -  -ip  u,oACd  /a(x)e 


<e 


|zlsl  +  e 


l*lpl)d* 


i  w  Aa 
2  o  o 


/k(-d.  +  ix  cosB) 

a(x)e  cos(kSD  sinfJ)dx 


(C.7) 


F5  "  1P  'W  “oACD  J  xa(x)e 


k(-dj  +  ix  co s 6)  ikSD  sin 8  t  -ikS^  sing 


(e 


n  # 

zls>  +  e  lzlpl)dx 


+  i 


w  Aa  U  f 
o  o  J 


k(-d,  +  ix  cos$) 

xa(x)e  cos(kSn  sin$)dx 


(C.8) 


where  the  Integrals  are  along  both  hulls. 
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f 


‘is 

zlp 

*10 


i3  -  x£5  -  (v(x,-SD(x)(-d1(x)) 
£3  -  x£5  -  Cv(x,SD(x),-d1(x)) 


Two  types  of  terms  are  present:  ones  which  include  CD,  a  cross  flow  drag 
coefficient,  and  ones  which  include  aQ,  a  body  lift  coefficient. 


CROSS  FLOW  DRAG 

In  investigating  the  force  F  acting  on  fixed  cylinders  in  oscillatory  flow, 
the  Morlson  equation  has  been  applied 

dU 

F  ■  C.M.  if  *  T  <c-9) 

is  an  inertia  term,  Aq  is  the  characteristic  body  area  normal  to  the  velocity 
vector,  U  is  the  fluid  velocity,  and  d  is  the  width  of  the  body  transverse  to  the 
flow.  Note  that  the  second  terms  in  Equation  (C.9)  and  Equation  (C.l)  are  equiva¬ 
lent.  In  an  experimental  investigation  of  two-dimensional  cylinders  and  flat 

23 

plates  in  sinusoidal  flow,  Keulegan  and  Carpenter  showed  that  C^  and  CQ  can  be 
represented  as  functions  of  the  parameter  U^T/d  where  is  the  amplitude  of  the 
harmonically  varying  velocity  and  T  is  the  period  of  oscillation.  This  parameter, 

K  ,  is  referred  to  as  the  Keulegan-Carpenter  Number.  Additional  experimental  work 

c  5 

on  circular  cylinders  has  been  reported  by  Sarpkaya.  He  introduced  an  additional 

2 

parameter,  6^, which  is  defined  as  d  /vT  where  V  is  the  kinematic  viscosity  of 
water. 

The  potential  for  applying  this  data  to  SWATH  ships  is  evident,  since  the 
lower  hulls  of  SWATH  configurations  are  typically  circular.  Experimental  data 
for  the  zero  speed  heave  added  mass  and  damping  coefficients  for  eight  configura¬ 
tions*  is  available.  Utilizing  one  circular  lower  hull,  three  strut  variations, 
and  two  draft  variations  for  one  strut  were  constructed.  Three  elliptical  sections 


^Reported  by  Stahl  in  an  NSRDC  Test  and  Evaluation  Report  of  limited  distribu¬ 
tion. 
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with  the  middle  strut  and  the  middle  draft  were  constructed.  The  results  from  the 
oscillation  experiments  showed  strong  dependence  of  damping  on  wave  amplitude  and 
some  dependence  of  added  mass  on  wave  amplitude. 

In  applying  Sarpkaya's  data  to  SWATH  ships,  results  for  CD  are  linearly 
Interpolated  as  a  function  of  both  Kc  and  8^.  In  addition,  for  each  8^»  results 
are  linearly  extrapolated  to  a  zero  value  of  K£.  For  large  Kc,  results  are 
assumed  to  be  equal  to  the  value  for  the  largest  K£.  However,  the  values  of  Kc 
encountered  for  the  experimental  conditions  are  small.  For  small  the  value  of 
is  2.0  which  corresponds  to  the  potential  flow  case.  Consequently,  Sarpkaya's 
data  do  not  alter  the  added  mass  results. 

In  Figure  11  typical  correlation  is  shown  for  the  case  where  Sarpkaya's  data 
are  applied  directly.  The  open  symbols  represent  experimental  results  and  the 
solid  symbols  represent  predicted  results.  A  different  symbol  is  used  for  each 
amplitude  of  oscillation.  Heave  damping  nondimens iona lized  by  mass  times  frequency 
is  plotted  versus  wave  number  nondlmens iona lized  by  2/(strut  beam).  Clearly,  the 
damping  is  over predicted. 

It  is  reasonable  that  a  circular  section  with  strut  will  have  less  viscous 
damping  than  a  circular  section  without  strut.  Based  on  the  correlation  for  the 
eight  configurations,  a  modification  factor  was  developed  to  alter  CD  for  applica¬ 
tion  to  SWATH  sections.  There  are  certain  limiting  conditions:  for  a  circular  sec¬ 
tion  without  strut,  CD  should  equal  Sarpkaya's  C^;  for  a  section  where  the  strut 
thickness  is  the  same  as  the  major  axis,  CD  should  be  close  to  zero.  An  expression 
which  satisfies  these  conditions  and  which  Includes  a  factor  to  alter  results  for 
ellipses  is 


C 


D 


JL 

a  a  / 


SARPKAYA 


(C.10) 


where  B  is  the  beam. 

Predicted  results  for  the  eight  configurations  in  Stahl's  investigation  are 
given  in  Figure  12.  For  elliptical  sections,  CB  is  evaluated  using  Sarpkaya's  data 
for  a  section  with  the  diameter  equal  to  the  major  axis.  The  predicted  potential 
flow  contribution  also  is  given  on  each  figure.  Although  the  experimental  and 
predicted  results  do  not  agree  precisely,  the  results  are  generally  encouraging. 
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BODY  LIFT 

In  Equations  (C.2)  through  (C.8)  terms  of  the  form 


-j-  aQU  J  x*a(x)dx,  i  -  0,  1,2 


appear.  These  terms  are  due  to  lift  on  the  body,  and  follow  from  Equation  (C.l). 
In  other  nomenclature  commonly  used  for  hydrodynamic  coefficients,  B,^,  B^^, 

and  B^  can  be  related  to 

-  UL2Z\  -  UL3M',  -  -E-  UL3Z\  and  -  UL4M' 

2  w*  2  w  2  q*  2  q 

*  *  *  * 

respectively.  That  is,  B^*  853*  835*  and  B,-^  can  be  equated  to  -Zw,  -Mw,  -Z^, 
and  -M^,  respectively.  In  Reference  24  it  is  suggested  that  submarine  data  be 
used  to  define  the  body  lift  expressions  for  SWATH.  This  approach  was  not  pursued 
by  Lee  but  is  pursued  here.  The  body  lift  components  in  Equations  (C.2)  through 
(C.6)  will  be  replaced  by  expressions  related  to  the  submarine  hydrodynamic 
coefficients  and  an  expression  for  aQ  will  be  developed  for  use  in  Equations  (C.7) 
and  (C.8). 

Dempsey  developed  semlempirical  expressions  for  vertical  plane  hydrodynamic 
coefficients  for  submarines.  For  a  bare  hull  submarine  the  expressions  are: 


Z^  *  2.439[(k2  -  k1)m,]4/3  -  m' 


(C.ll) 


*  0.793(k2  -  kj)m' 


(C. 12) 


Z'  -  -0.207(1^  -  kj)m' 


M’  -  (Z')2/Z' 

q  q  v 


(C. 13) 
(C.14) 


where  m'  ■  2V/L3  with  V  equal  to  the  bare  hull  volume  and  L  equal  to  the  ship 
length.  Lamb's  hydrodynamic  coefficients,  kj  and  kj,  are  given  by  : 


*Given  in  a  DTNSRDC  report  of  higher  classification. 
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where  a  is  the  maximum  of  the  horizontal  axis  of  the  lower  body, 
max 

In  applying  these  equations  to  SWATH  configurations,  the  obvious  approach  is 
to  evaluate  the  submarine  expressions  for  one  SWATH  lower  hull  and  double  the 
value.  For  the  SWATH  6A,  comparison  between  predicted  and  experimental  results 
showed  that  the  calculated  results  were  lower  than  the  experimental  results.  Con¬ 
sequently,  in  evaluating  Dempsey's  expressions  for  SWATH,  the  volume  is  taken  to 

Include  the  volume  of  the  strut  as  well  as  that  of  the  submarine-like  lower  hull. 

* 

Dempsey's  expressions  could  not  be  applied  directly  for  B^*  When  it  was 

taken  to  be  -2Mw,  where  the  factor  of  2  is  used  to  include  the  effect  of  both 

hulls,  it  differed  in  sign  and  significantly  in  magnitude  from  the  results  for 

the  experiment  documented  in  Reference  12.  However,  the  difference  between  the 

experimental  results  and  the  sum  of  the  potential  flow  and  cross  flow  drag  com- 

* 

ponents  is  approximately  equal  to  -2Z^.  Therefore,  is  taken  to  be  equal  to 
-22^.  That  is,  is  taken  to  be  equal  to  Bjj.  This  is  consistent  with  Equation 
(C.3)  in  the  theoretical  development. 

Appropriate  body  lift  components  for  and  must  be  determined.  Refer¬ 
ence  to  Equations  (C.2),  (C.3),  (C.S),  and  (C. 6)  Indicates  that  B^  and  as 
well  as  8^2  end  are  similar  In  form.  Analysis  of  data  taken  for  the  SWATH  6A 
during  the  experiment  documented  in  Reference  12  Indicates  that  for  that 
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A  ^ 

configuration,  at  leant,  C35  can  be  approximated  by  -2UZw  and  that  can  be 


approximated  by  -0.278  UMv< 

In  order  to  evaluate  the  excitation  components,  some  value  for  a  must  be 

o 

chosen.  Several  possible  relationships  follow  from  Equations  (C.2)  through  (C.4) 
and  the  corresponding  Equations  (C.ll),  (C.13),  and  (C.14).  For  the  SWATH  6A, 
Equations  (C.2)  and  (C.ll)  yield 


L2  Z ' 
_ w 

/  a(x)dx 


0.068 


(C. 15) 


Equations  (C.3)  and  (C.13)  yield 


L2  Z ' 


/x  a  (x)dx 


2 —  =  3.98 


(C. 16) 


Equations  (C.4)  and  (C.14)  yield 


L3  M' 


/x2a(x)dx 


0.127 


(C. 17) 


The  value  which  is  given  in  Equation  (C.15)  is  close  to  the  value  of  0.07 

2 

suggested  in  Reference  22  and  used  by  Lee  .  In  Reference  22  the  experimental 
data  which  produced  this  value  was  for  the  heave  force.  Thus  it  is  reasonable 
that  it  should  agree  with  the  value  from  Equation  (C.15).  The  other  values  differ 
substantially  and  are  a  measure  of  the  difference  between  using  Equations  (C.2) 
through  (C.8)  with  a  constant  value  of  aQ  to  represent  the  body  lift,  and  using 
the  above  relationships  derived  from  the  submarine  hydrodynamic  coefficients  to 
represent  the  body  lift  component. 

These  three  values  for  aQ  were  substituted  in  Equations  (C.7)  and  (C.8). 

Correlation  with  experimental  results  for  heave  exciting  force  and  pitch  exciting 

12 

moment  for  the  6D  and  the  SSP  KAIMALINO  Indicates  that  the  best  correlation 
results  from  using  Equation  (C.15)  to  evaluate  aQ. 

FIN  CONTRIBUTIONS 

Contributions  to  the  added  mass  due  to  the  presence  of  stabilizing  fins  are 
given  in  Lee's  development.  In  addition,  a  development  similar  to  that  for  the 
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body  results  In  fin  contributions  due  to  cross  flow  drag  and  lift  in  the  damping 

coefficients.  The  resulting  equations  require  evaluation  of  a  cross  flow  drag 

coefficient,  Cj^,  and  the  lift  curve  slope,  ,  for  the  ith  fin.  Lee  suggests 

using  1.2  for  CQ.  and  uses  a  technique  for  evaluating  C.  which  Includes  the  effect 

25  1X11 

of  the  body  on  the  fin  and  the  fin  on  the  body.  Alternate  approaches  for 
evaluating  these  coefficients  are  possible. 

Data  in  Reference  6  are  used  to  evaluate  These  data  are  given  as  a  func¬ 

tion  of  Kc,  the  Keulegan-Carpenter  number,  UmT/d.  They  are  presented  in  Figure  10, 
in  which  values  for  a  flat  plate  range  from  2.5  to  6.0.  Kc  is  evaluated  using 
a  value  of  twice  the  span  for  d  with  as  the  amplitude  of  the  harmonically  vary¬ 
ing  velocity.  T  is  the  encounter  period  as  for  the  values  of  the  body. 

As  noted  by  Dempsey, ^  work  by  Whicker  and  Fehlner^  showed  that  for  low  aspect 
ratio  control  surfaces  in  the  free  stream,  experimental  results  agreed  with  this 
expression 

c  =  _ 1 . 8tt  AR 

10  1.8  +  (AR2  +  4.0)1/2 


where  CT^  is  the  slope  of  the  lift  coefficient  with  respect  to  angle  of  attack  a 
in  radians  and  a  rectangular  planform  has  been  assumed.  AR  is  the  effective 
aspect  ratio  which  for  the  SWATH  case  is  given  by 


c 


where  s  is  the  span  and  c  is  the  chord  of  the  lifting  surface. 

A  source  of  Interference  on  lift  is  the  effect  of  the  hull  on  the  stern  foil. 
Dempsey1^  developed  a  semlempirical  expression  for  a  correction  to  for  stern- 
planes 
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0.2556 
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(a  + 


//  a  +  2s  V 

-  0.1612^ 

W  a  / 
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1/2 


-  0.6366  sin 
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(a  + 


0.4015 
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Application  of  her  expression  to  the  SWATH  6A  forced  osclllat  on  results  in 
Reference  12  for  one  and  two  fins  yielded  poor  correlation.  Analysis  of  the 
respective  geometries  suggests  why  this  correlation  generally  will  not  be  appro¬ 
priate  for  SWATH  hulls.  The  stemplanes  of  a  submarine  are  very  near  the  stern  so 
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that  the  wake  of  the  submarine  will  mask  the  fin  and  reduce  lift.  For  SWATH 
configurations  with  the  aft  fins  set  near  the  parallel  mid-body  the  effect  of 
hull  outreach  will  be  insignificant.  The  correction  is  included  when  the  ratio 
of  the  maximum  transverse  axis  to  the  transverse  axis  at  the  quarter  chord  is  at 
least  2. 

Fin-fin  interference  also  is  a  source  of  change  in  lift.  Experimental  work 

g 

by  Lloyd  investigated  the  frequency  dependent  interference  effects  on  lift  due 

to  the  presence  of  more  than  one  appendage.  In  this  experiment  two  identical 

fins  with  an  aspect  ratio  of  2.0  were  attached  to  a  board.  For  separations  of 

10,  20,  and  30  times  the  span,  the  fins  were  oscillated  at  several  frequencies 

q 

at  amplitudes  of  oscillation  of  10  to  20  degrees.  Cox  and  Lloyd  corrected  this 
data  for  boundary  layer  effects  and  tabulated  the  data  in  terms  of  the  ratio  of 
the  lift  on  the  aft  fin  to  the  lift  on  the  forward  fin.  In  tabular  data,  chord 
was  used  as  a  nondlmensionalization  factor.  For  application  here,  the  outreach 
is  of  Importance  and  consequently  parameters  given  in  Reference  9  have  been  re- 
nondlmensionalized  using  span.  Reference  to  Table  4  shows  that  over  the  range  of 
values  tested  the  effect  of  interference  ranged  from  a  59%  reduction  to  a  15% 
increase  in  comparison  with  the  lift  of  a  single  fin. 

These  experimental  data  apply  to  the  case  where  both  fins  are  the  same  size. 
This  is  not  the  usual  case  for  SWATH  configurations.  In  the  absence  of  data  for 
various  fin  size  ratios,  the  interference  factor  associated  with  the  span  of  the 
forward  fin  was  used  to  modify  the  free-stream  lift  curve  slope  of  the  aft  fin. 
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Figure  11  -  Experiment -Fred let  ion  Comparison  for  Damping 
Coefficient  with  No  Correction  for  C^ 

^Reported  by  Stahl  in  a  report  of  limited  distribution. 
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Figure  12  -  Experiment -Prediction  Comparison  for  Damping 
Coefficient  for  Two-Dimensional  Cylinders 


Figure  12  (Continued) 


Reported  by  Stahl  in  a  report  of  limited  distribution. 
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Figure  12  (Continued) 
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Figure  12d  -  Coefficients  for  Configurations  with 
Elliptical  Lower  Hulls 


Reported  by  Stahl  in  a  report  of  limited  distribution. 
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